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ANNALS OF MATHEMATICS. 



Vol. III. October, 1887. No. 5. 

ON THE CONDITIONED COOLING AND THE CUBICAL CONTRACTION 

OF A HOMOGENEOUS SPHERE.* 

By Mr. R. S. Woodward, Washington, D. C. 

13. In the former paper (Annals of Mathematics, Vol. III. No. 3) we have 
di.scussed the problem of the "free cooling" of a homogeneous sphere. In the 
present paper it is proposed to consider what we Have, for brevity, called the case 
of "conditioned cooling of a homogeneous sphere," or the cooling of such a 
sphere in a diathermous medium. In addition, we shall give an application of 
the theory of cooling in determining under certain limitations the consequent 
cubical contraction in large spheres, this problem being of special interest with 
reference to the secular contraction of the earth. 

The physical circumstances, and the mathematical features of the solution 
of the problem of conditioned cooling have been very fully and ably set forth by 
Fourier.t Poisson,| and Riemann.§ So far as the purely analytical requirements 
are concerned, their works, e.specially that of Poisson, will be found exhaustive. 
We need, therefore, only to refer to them for a complete exposition of the ele- 
mentary considerations which lead to their final result, our object being to adapt 
that result to the requirements of the computer, having in mind always an ap- 
plication of the theory to the earth. It will be desirable, however, for the sake 
of clearness and definiteness to make a formal statement of the fundamental as- 
sumptions and relations involved in the steps of the investigation which precede 
those we have in view. 

* Since this paper is es.sentially a continuation of the one on " The Free Cooling of a Sphere," pub- 
lished in Vol. III. No. 3 of the Annals of Mathematics, the section numbers and those designating the 
principal equations are made continuous with the corresponding numbers of the earlier paper. The nota- 
tion introduced in the latter has been adhered to as far as possible. 

fThiorie Analytique de la Chaleur. Paris, 1822. Analytical Theory of Heat, by Joseph Fourier. 
Translated, with notes, by Alexander Freeman. Cambridge, 1878. 

|Th4orie Math^matique de la Chaleur. Paris, 1835. 

J Partielle Differentialgleichungen und deren Anwendung auf physikalische Fragen. Vorlesungen 
von Bernhard Riemann, fiir den Druck bearbeitet und herausgegeben von Karl Hattendorff. Dritte Auf- 
lage. Braunschweig, 1882. 
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14. The assumptions made are four, viz. : first, a constant coefficient of dif- 
fusion for the entire mass of the sphere ; second, an initial uniform excess in tem- 
perature above that of the surrounding medium of the entire mass of the sphere; 
third, a rate of flow of heat from the surface of the sphere to the medium pro- 
portional to their difference in temperature ; and fourth, a constant temperature 
for the surrounding medium. 

As in the previous paper, let r be the radius of any spherical shell whose 
thickness is dr, u the excess in temperature over that of the surrounding medium 
of this shell at the time t after the initial epoch, and (f the coefficient of diffiision. 
Then these quantities must satisfy the partial differential equation 

and the solution of the problem consists in determining the integral of this equa- 
tion, {ru) being regarded as the unknown quantity, subject to certain relations 
which we proceed to set down. If «„ be the initial uniform excess of the tem- 
perature of the sphere over that of the surrounding medium, we must have 

ru = rUf, for t = o. (30) 

In order that u may not be infinite when r = 0, we must have 

rw = o for >- := o. (31) 

It remains to express the characteristic condition of the problem, namely, 
that which defines the emission of heat at the surface of the sphere. If k be 
the coefficient of conductivity for the sphere, the quantity of heat which passes 
across a unit of its surface in an element of time is expressed by the value of 

— k-fr- when r = r„, the radius of the surface of the sphere. On the other hand, 

this quantity of heat is equal to that which the surrounding medium can carry 
away, or according to our third and fourth assumptions, to Hu, H* being a con- 
stant and u the temperature corresponding to r = r,,. Hence we have 



or, if we put H\ k = h, 



But since ; 



(rii) 





-4"=*. 


J- -]- ku = o for r = 


du 


I d (ru) u 


3r~ 


r dr r 



* This quantity is called by Sir W. Thomson the emissivity . It is the quantity of heat emitted per 
unit of time, per unit of surface, per unit of excess of temperature ; and may be due to conduction, con- 
vection, or to both. See \']X, Article Heat, Encyclopedia Britannica. 
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which, substituted in the above relation, gives 

). ' + (rk — i)u = o for r = r^. (32) 

Now a particular integral of (29) satisfying the relation (31) is 

ru = C'e-''"''sinXr, (33) 

in which C is a constant, and ). a factor such that 

^^^ =C'Xe-'''"''cosXr 

dr 

and u= C'r-^e-''"" smkr 

will satisfy (32). Substituting these values in (32) we have X defined by the fol- 
lowing transcendental equation : — 

r^^X cos r^X -^ (rf,h — i) s\n r^X = O. (34) 

This equation has an infinite number of unequal real roots,* and since any one 
of them will satisfy (33), the complete integral of (29) must consist of a series of 
terms of the form (33). That is, 



« ^ 00 



2> 



ru= I C'^e-"''^"' s\nX„r. (35) 

» = I 

In this, Xi, X^, X^, . . . are the successive roots of (34). The constants C\, C'^ 
C\, . . . must be such that (35) will satisfy the initial relation (30); i. e. making 
tr=o in (35), the relation 

ru^ = I C'n sin X„r. 

The process of determining C„ from this equation consists in multiplying both 
sides by sin X „rd r gind integrating between the limits o and r,,, account being 
taken of the relation (34). We shall not go through the steps of this process, 
referring the reader to Poisson and Riemann for the details, but shall pass imme- 
diately to one of the usual final forms for (35). This form isf 
« ^ 00 

«=i (36) 

15. The expression (36) satisfies all of the relations (29) to (32) , because each 
of its terms satisfies them, and is the complete solution of the problem. But, 
though it meets all the demands of pure analysis, its form is quite repulsive to 

* See any of the works 6rst cited above. 

f See Riemann, Differentialgleichungen, J69. Other equivalent forms are given by Fourier, Poisson, 
and Riemann. 
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the computer or to any one who may desire to use it, considering that the quan- 
tities ^,„ must be determined from the transcendental equation (34). Some at- 
tempts have been made to simplify this expression, but the results are not satisfac- 
tory, inasmuch as they require either a time so great that we may neglect all terms 
of the series after the first, or a very small radius of the sphere, or a very great 
radius of the sphere. No attempt, so far as we know, has been made to trans- 
form the series into another more converging or tractable, though Poisson has 
summed it to a limited degree of approximation, as we shall presently explain. 
Neither is the improvement we have to suggest entirely satisfactory, for it is re- 
stricted to the case where r^h of (34) is greater than unity,* but this is the case 
with which we are most particularly concerned and that to which we now direct 
attention. 

16. The quantity h in (34) is of — i dimension in length, f and hence rjt is 
of zero dimension or an abstract number. 

For brevity, let us put 

x=\\ {rjt) and j/„ = r^^X,,. 
Then equation (34) gives 

t^ny^ = --^. (37) 

Likewise, if we make the same substitutions in (36), it may be put in the more 
compact form 

ru-zru I— (i-^)'+^V sity^ ^-»"*..'sin; r ix'Ss 
»^ I 
Now it will be observed that when x^=o,Vn^^ nn, X„ = nn j r„, and 



hence (38) becomes 



ru 



This agrees with equation (10), which expresses the solution in the case of 
free cooling. The case of conditioned cooling, therefore, approximates to that 
of free cooling in proportion to the smallness of the number x. This number 
may be small by reason of the largeness of either r„ or h, or of both. From cer- 

* It is obvious, however, that a slight modification of the process followed in the text may be used in 
treating the case where r„ A <] i. 

fSee Fourier's Section IX. chap. II. p. 130, Freeman's translation. Our h is the same as Fourier's 
hi Xof the page just cited, but the same as his h defined on p. 268. 
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tain observations of underground temperatures in the garden of the Paris Obser- 
vatory, Poisson concluded that for the earth in the vicinity named, ^ is a little 
greater than i, using the metre as space unit* This would make x= i \(rji) 
about I / 6,300,000. Whether this value is nearly correct as applied to the whole 
surface of the earth may be doubted, but x might be 10,000 times as great with- 
out materially affecting the applicability to the earth of the process we shall 
develop. The correctness of that process only requires that x <^\. 

It is proposed, now, to expand the second member of (38) by Maclaurin's 
theorem in a series of ascending powers of ,r. As a preliminary, though not es- 
sential part of this work, we may expand y^ of (37) as a function of x and arrive 
thus at a formula of use in a future numerical application. 

17. When X = 0, (37) gives y^ = nn, and this is the first term in the expan- 

d'y 
sionof y„. To get the coefficients of the terms involving powers of x, we need -j^ 

for x ^o,i being any positive integer. If we dififerentiate (37) once, and then 
eliminate sec^^„, the result is 

If now we put Y=y^-j^ , and dififerentiate the above equation successively, we 
find that for jr =: o 

^" = (^-2)^^"-(^- 1)0-2) ^j^. (39) 

By a duplicate application of the theorem of Leibnitz, we can derive a gene- 
ral expression for the (/ — 3)th diflferential coeflficient of y with respect to x, but 
its usefulness hardly justifies the space required by its rather complex character. 
We shall simply write down the first six values of (39). They are 



2;- " 
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i= I, 


= 




= 2, 


= + 2{nnf 




= 3, 


= — 8 {HTzf 




= 4, 


= — 24 (mzf 




= 5, 


= + 384M' 




= 6. 


Hence by Maclaurin's theorem we have 






y„ = n7:{i—x^ ^n^nV — 


i«= 


't^x' — . 


»Thtorie de la Chaleur, p. 502. The value given is 1 


i = . 


h = I.05719. 



)• (40) 
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Neglecting terms involving powers of x higher than the 4th, this equation may 
be written 

(41) 
= nn — nnx { i — ^«V;r* ( i — x)}, 

one of which forms will be found convenient for logarithmic computation of the 
smaller values of j„. For the larger values, the series (40) will, in general, prove 
too slowly converging, but it is not necessary to calculate such values. Our only 
use for the smaller values is a numerical verification. 

18. To develop (38) we observe that each of the seven factors under the sign 
of summation, is a function of ;t;. We might, therefore, expand these factors sep- 
arately and then find their product. On the other hand, we might get the differ- 
ential coefficients of (ru) with respect to x by direct dififerentiation. We shall 
adopt a combination of these two processes. 

First, develop the product x~^ sin_j/„. By means of (40), we readily find 

x~^ sin j„ = — nn cos njz (i — ^nh:V + \t^7j^x^ + . . . ). 

Then, if we put y„ = tiny, where y is the series within the parentheses of (40), we 
see that 

j/„-^;r-'sinj/„ = ^~ '^"^' .y~'(l — \n^T^x^ -\- \n^i^:^ + . . .). 

Now, for brevity, let 

X^= l-x +x'y^\ 
and X^ = {i— ^nVx^ + ^k'tt'x^ + . . .). 

Then, if we denote the «'* term of the series under the sign of summation in (38) 
by T„, and dififerentiate log T„, the result is 

T^^-~^y Tx^^' -d^-^^ ^+^' ^ 

-..^..„§ + .cot;„.§. 

But from (40), we have when x^o,y= i , and -r- = — i . Also, since A ^ -^ 

dx r^ '■ 

. nn d^„ MTT 

" ~ r„ dx~ r„ ■ 



Therefore, when ;r = o, we find 



'^• = ^- [■ + -[■ 

The succeeding complete values of the differential coefficients grow rapidly in 



^_7;|_,+2^i^— J --nncotnn-J, (42) 
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complexity. Their derivation, however, i.s tedious rather than difficult. It will 
suffice, for our purposes, to write down the special value of the second differen- 
tial coefficient, viz. : 



f 4^«v[,+:-;]+34^]-l 



9a^ 



= T„ 



+ ^■[^1 



i'^"" + "»'"• (7^]"] 



cot HTT 



(43) 



Now we observe that when x = o, 



T = 



71 — tm cot nn — = 
fo ''0 



Cann 






dT * 
r —^ 

dr ' 



e -"" ^"^ I '0)' sin nn- 



2t 



dt ' 






4r.*4=] 



cot ««■ — =+ 4^/ 



9r9t' 



(^) 
(^) 

(^) 
(/) 



For brevity, put 5 for the value of the second member of (38) when x = o, this 
value being the same as that given by (10), (19), or (20). Then, observing the 
effect of the sign of summation in (38), equations (42) and (43) in connection 
with the equivalents in (a) to (/) give 



9(ru) 

~W"' 



3\ru) 
3x^ 






4S — 2r 



3S 

9r ' 



9S 






3r^ 



3^S 



(44) 



(45) 



*To get the signs of the second members of (J>) and (/), one must note that STnl Qr and dknj dx 
have opposite signs. 
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Finally, by Maclaurin's series we get 



Sx 



dx^ 



(46) 



19. The partial differential coefficients involved in the second members of 
(44) and (45) are easily calculated. If we write, for brevity, 



Q' = I I 



« =00 i 
K= I J 






(2 Wo« — w„ — m) e - <='""'" -"•«-"•)'' 
(2»«„« + Wo + w) f-C*".- + «» + '»)'' 
(2 W9« — m^ — ntfe- <^'"»" - "» - "■>' 

(2W9« + Wo + '«)'' ^ ~ *^°" "^ '"° "^ ""' 
(2 Wo« — m^ — mfe- <^'"»'' - "» - "*>' 
(2Wo« + w„ + w)' ^-C2»„» + », + ».). 



Wo and m being defined by equations (18), we get from (23), (24), and (25) 



(47) 
(48) 
(49) 

(SO) 

(51) 
(52) 
(53) 
(54) 
(55) 



These will be applicable for all but the latest stages of cooling, and the corres- 
ponding derivatives from (10) will become applicable when these fail to be suffi- 
ciently converging. 

The solution for conditioned cooling is thus made to depend on that for the 
ideal case of free cooling ; and the expression (38), which is practically useless 
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(on account of its slow convergence) for small values of the time is transformed 
into (46), which is available for any value of the time. We have limited the ex- 
pansion to terms of the second order in x, believing that this is more than suffi- 
cient for any application of the theory to the earth; but the higher terms of the 
series may be easily derived if needed. 

20. To complete the solution, it remains to show how to determine the tem- 
perature at the centre of the sphere. When we make r = o, the expression for 

u from (46) assumes the ambiguous form - ; but the usual process of differentia- 
tion gives 



From equation (26) we have 



« ^ 00 



"-»']. (57) 



The differential of (44) gives 



From (55) we have 



= 2 ^- — 2/-^ (58) 

or 3r 3rdt 



^'^. = .-%7)(S'-^e"'), (59) 



in which, of course, r must be made zero for the special value in (58). 

The special value of the differential of (45) with respect to r, required by 
the second term of (56), may be easily written down and worked out numerically 
for any assumed data, but the paucity of practical applications makes this work 
unneces.sary.* 

*In note Cof the supplement to his great work (see full title, p. 129) Poisson has summed the series 
(38) to terms of the first order in x. He does not separate his result into the term independent of and 
that dependent on x, but gives an expression for the sum of the two terms. He deduces several special 
results dependent on the assumption (a) that jr = o, or (/>) that the strata whose temperatures are consid- 
ered are near the surface of the sphere. Among these results is an expression for the temperature at the 
centre of the sphere under the assumption that x = 0. This expression agrees with our expression (26) or 
(57) ; but the casual reader would not observe the agreement; and, curiously enough, Poisson himself was 
unaware that the difference of two definite integrals in his expression is exactly zero. Using our notation,. 
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21. Since no set of formulae presenting any degree of complexity, or requir- 
ing many steps in their derivation, can be regarded as entirely trustworthy until 
they have been subjected to the test of actual application, we propose in the pres- 
ent and next sections to work out in detail a numerical example. For this pur- 
pose we shall assume such a value of the time that we can compute ru from (38) 
as well as from (46). The equivalence of the two forms will thus be tested, as 
well as the correctness of the signs of the several terms in the second members 
of(44)and(45). 

We shall assume m^ = r^j 2ai/t = tt/ 2, a value used also in the illustrative 
example of § 1 1 . This gives 



Ian J ' 



which, for the earth, if we use Sir W. Thomson's value of a^, is, as stated in §3, 
about 100,000,000,000 years, With this value of the time, the first three terms 
of (38) will suffice to give ru in terms of r^^u^ to seven places of decimals. We 
shall also assume x= 10 ~' and r = ^r^ and compute first by means of (38). 

22. To begin with, we need the first three roots of the transcendental equa- 
tion (37). These we get from the first of (41), which gives with sufficient accu- 
racy log^„ = log {nn) -\- log (i — x) -\- \iJL {nnf x^; /J- being the modulus of the 
common system of logarithms. This gives 

log:>'i = 0.4967154, 

log J2 = 0.797745, 

logj»'j = 0.973837. 

For use in sin_y„ and sin X^r, it is be.st to compute the supplements {nz — _>/„). 

he has for the fall of temperature at the centre of the sphere (see p. 52 of the supplement) 

CO « = « n =- « + jrto 

y Kid n^l n.— 1 — »rto 

while our equation (26) gives only the middle term of this expression. But if in the last definite integral 
we replace znii^n -|- 2 by 0, the result is 

» =— 00 flio (2ii + 1) 3nio 5»«o 

2' / e-"de= / e-'"'dz-\- / e-"dt:-\- . . . 

» = 1 mo (2n — 1) mo Sm© 

00 

= / e-"dz, 

which is identical with the first definite integral. Hence the above expression is the same as (26). Pois- 
son, after reaching his equation, observes that 

+ '"0 

— mo 
and proceeds to draw some conclusions which are qualitatively correct, but not susceptible of definite nu- 
merical evaluation. 
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These may be accurately computed from the second of (41), which gives in sec- 
onds of arc 

;t— J, = 64;'.'97;87, 

2ff— j/2= 1295.983. 
?>'^—J'i= 1943- 94- 

These and the above values of log/n will be found to satisfy (37). 

We can now compute all the factors in (38). The logarithms of the second 
algebraic factor and those of the factor x~^ sinj„ are found to be 

log XiX^ - ' log ;ir- ' sin /„ 

9.9995655 — ID, 0.4971479, for « = I, 

9.999565 — 10, 0.798170, = 2, 

9.999566 — 10, 0.974252, = 3. 
Since ^„ =zy^lr^, we have with our a.s.sumed value of/ 

g-n'l\„' __ ^-(.v„;t)» 

and hence log e-"'"'"' = 9.5665736 — 10, for n^= 1, 

= 8.266294 — 10, ^ 2, 

^6.099162 — 10, ^3. 

Likewi.se, for the last factor we have sin X„r = sin_y„ rj r^ = sin ^y„. Hence 

log sin ^j/„ = 9.9999995 — 10, for «= I, 

= 7.798170 — ID, = 2, 

= 9-99999S» — 10, = 3- 

Summing the logarithms of the factors of the several terms we get 

log No. 

I. St term, 9.0698557 — 10, +0.1 1745073, 
2d 4.965681 — 10, + . 0924, 

3d 5 125300,, — 10, — . 1334. 

The sum of the terms is o. 11 744663. Multiplying this by 4 we get as the result 
from (38), 

« = o.4697865«„. 

23. The details of the computation by equation (46) are as follows : — 
The value of twice the first term, S, has already been computed in the ex- 
ample of §11 ; i. e. 

5 = o.234i732r„«„, 
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using the larger of the two values there given, since a more extended computa- 
tion shows that to be the more precise value. 

We next need the values of the Q's which enter the second members of (51) 
to (ss). Since hy {18), m ^ r 1 2a\/t, and since for this example we have assumed 
m^^=r^l 2a]/ 1 ^ tt / 2 and r = ^r^, we have for use in (47) to (50) 

m^ — m^ \iz, fftf^ + M = f n-, 
3W„ — m = ^7r, 3w«„ + tn= |7r, 
5w„ — m ^^7z, 5wz„ + »« = ij-K, 



Hence we have, to terms of an order which may be neglected, the following sys- 
tem of equations in which the exponential quantities are factors respectively of 
the several coefficients in the vertical rows : 

Q' = I ^-OT=+ I ^-<i'^)' + I ^-(i'^)'+. . . 
Q" = In - |;r + f;r - . . . 

Q"'=iiny +{inf +{iny +. . . 

&-^i\ny -iinf +{i7:f -. . . . 

These give, by an easy computation, 

Q' = + 0.54352, 
Q" = + 0.41469, 
Q'" = + 0.35443, 
giv ^ + 0.21069. 

Introducing these factors and the values of i and r in equations (51) to (55), we 
find the following equivalents, in which r^u^ is a common factor of the second 
members : — * 

r^= + 0,01832 r„2<„, 
2/ -^ = — 0.46793, 

4f'^ = + 0.928. 

♦The numerical values of all these differential coefficients have been checked by means of the equiv- 
alent forms obtained from equation (10) . 
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By means of the last values and the above value of 5 we get from (44) 
Likewise from (45) 



9 (ru) 

-^= + 0.72042 r„?<„. 



-g^=— 0.883 r„«„. 

With these values and with x = 10 ~', we have from (46) 

u = 2«„ (0.2341732 + 0.00072042 — 0.00000044) 

= 0.4697864?/,,. 

This agrees with the value obtained from (38) as closely as 7-place tables of log- 
arithms will permit, the discrepancy being one unit in the last place. 

It is scarcely necessary to remark in this connection that the precision of the 
preceding computations is of little importance except as a test of the correct- 
ness of the formula; derived and as an illustration of the process which might 
be followed if our data for practical applications were sufficiently trustworthy. 
For all applications to the earth, computation of results by the preceding for- 
mulae may be properly limited, in general, to three or four places of decimals, 
since defects in the best data at present attainable and probable inaccuracies in 
our fundamental assumptions (of §14), render those results uncertain even in the 
first place of decimals. For these reasons, also, and in view of the probability 
that the quantity x- is less than i / 1000, it is believed that the solution for free 
cooling is an ample approximation in all applications to the earth. 

24. In this and the following sections it is proposed to determine the amount 
of cubical contraction of a homogeneous sphere on the assumption of free cool- 
ing, and on the additional assumption that the coefficient of contraction is con- 
stant. These assumptions will give correctly the first approximation (or the term 
independent of jt) in the case of a homogeneous sphere of moderately high initial 
temperature, and probably as good an approximation as is at present attainable 
for the secular contraction of the earth.* 

If s denote the coefficient of cubical contraction and J V the total change in 
volume during the time t after the initial epoch, we shall have 



JV^ 47re I («„ — u) r^dr, (60) 



in which the symbols under the sign of integration have the same meaning as 

*The consideration of the mechanical effects of secular contraction forms no part of our present 
work, though one of the principal objects of the analysis of this and the preceding paper is to furnish a 
basis for the investigation of those effects. 
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heretofore. The fall in temperature, «„ — «, is given by equation (21) or (22), 
but for our present purposes it will be convenient to use another expression easily 
derived from (10). If in the latter, we make t first zero, and then t, and take the 
difference between the two expressions the result is 

r(«„ — «) = ^^^» I ^~ '^'"^' (i— g-^'t"" "■'-'') sin ««■-■ (61) 



TT n r^ 







Introducing this value of the fall in temperature of the shell considered, in (60) 
we get 

dV= ■ir^u^i I ^~ '^"^' (i — ^-«"C"r/r,)») r^^^ sin „^ ^ 
ft fc/ ^0 

» = I o 

But the integral in this equation is equal to 

r^ cos nn r/ ( — i ) " + ' 

n n t: n ' 

Therefore we have as the exact solution, under the assumed conditions, 

« = 00 
jy=^rlM Z \{i _^-»"(»T/-.)'),. (62) 

« = I 

When ( = 0, this expression gives J V= O, as it should. When / = 00 , the 
quantity under the sign of summation becomes 

and hence JV=^nr\e, 

which is also obviously correct. 

25. For large values of the time the series in (62) is evidently rapidly con- 
verging, but for small values the convergence is too slow for practical applica- 
tions. It is desirable, therefore, to transform (62) and derive an expression from 
which we can compute the contraction when i is small, or during the earlier stages 
of cooling. To this end we proceed, 

Referring again to the Gamma function of §4, and making the same substi- 
tutions for a andjS as in §5, the series in (62) becomes 
00 

•^ff* 7f—\ f^ ~ '''""''^" dx (cos nx -{-\ cos 2nx + ^ cos 3;r;i: + . . . ) , 

o 
in which x is merely a subject of integration and is, of course, not the x of 
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§§i6 to 23. Now, for brevity, denote the trigonometric series in the last expres- 
sion by (Tj, and let 

da^ d^a^ da^ 

dx " d:i^ dx 

From the series in §6 we see that 

a = l-!r ^-^ ^-^ — ; — 5, when f = i . 

"* I — 2g cos nx ^ g' " 

Multiplying by «^;tr, integrating between the limits O and x, and making ^= i, we 

get 

(T^ = ^T^ (x — i) for o < ;i; < 2. 

Multiplying the last equation hydx and integrating again between O and x, there 
results 

This gives the value of ff^ for all values of x from o to 2 both inclusive, or through 
the first period. For succeeding periods we have 

ffj = a^ [^ + ^ (;tr — 2f — ^{x — 2)] , for X between 2 and 4, 

= '^ [i + i (-^ — 4)^ — ^ (-^ — 4)] • *<^'' ^ between 4 and 6, 

Substituting these values of a^ ^^ '^he above integral, and writing, for brevity, 
m„ = r^fl za^t, as in equation (18), we get 
'2 

Je-"'"dx{^ + \x'-i^x) 
o 
4 
+ p-»'"'- dxli + i(x- 2f -\(X- 2)] 

2 

6 
+ J. - '"»-' dx\_\ + \(x- 4)^ -\{x- 4)] 

4 
^+ ^ 

Collecting and arranging the terms of this expression and performing the inte- 
grations so far as possible, it reduces to 

» = 00 2n ~\- 2 

— (— —} -h ff5 2' \ m-^ e-^^""'^' — 2mM (n -\- i) Ce"^'''' dx] . 



^r* — zm^iri ■< 



Substituting this value for the quantity under the sign of summation in (62) and 
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replacing x by s\m^, we get finally 

I I 



JF=8;i-A->„£ { 



ff/f, {2M -\- 2) 



^•(63) 



26. The series in the right-hand member of (63) converge with extreme sud- 
denness for small values of the time. For the case of the earth they converge 
so rapidly that their first terms alone will suffice for the most accurate computa- 
tions (one part in a million, say) when the time is as great as r//4^^, which, with 
Sir W. Thomson's value of a^, is about 273,000,000,000 years. Using this value 
of the time we have computed the amount of the contraction by both (62) and 
(63). The first three terms of (62) will suffice in this case for the precision at- 
tainable with 7-pIace logarithms. The second term in the first series of (63) or 
e~^^ j \/n amounts to but six units in the eighth place of decimals; and the sec- 
ond term of the second series, or 

4 
amounts to but twenty-eight units in the eighth place. Omitting other details of 
the computation, the results are 

JV=0.1S 80729 . 8n r„\e from (62) , 
= 0. 1 580728 . Snr„\e (63) . 

These results serve to verify the numerical equivalence of (62) and (63), and show 
that they, like (46), or (10) and the group (19) to (22), bridge completely the 
time interval of cooling and enable the computer to give a definite answer to any 
question falling within the scope of the assumed conditions. 

It is worthy of remark that the contraction just computed is about 

z(F=0.95.|;rr„X«; 

that is, in the time t^r^\^, a sphere will, under the assumed circumstances, 
attain ninety five per cent, of its possible contraction. The temperature at the 
centre of the sphere at the end of this time is (see §12) about one-sixth its initial 
value. 

In all practical applications to the earth, for any but the late.st stages of cool- 
ing, we may neglect the series in (63), which becomes then simply 



.rF=8..,v[f^^-g]. (64) 



From this it appears that in the earlier stages of cooling the contraction is pro- 
portional to the square root of the time. 



